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A . ovKRVin; 

Cur rojcarch durinq this tine |criod January I to J"ne ] , 1075 will be 
su.Tirarized undi'r the follovinq topicii (nee Section D) 

1) Systerr.n Feliability Optin; tion 

2) Failure Det« 2 Ction iMqoritb.ins 

3) Analysis of Ilon' ineir Filters 

4) Cosicjn of Conper.sator. s In.ror»-oratii;q Tine Delays 

5) riqital vTotr.peasator Cosiipi 

6) F'timition for Syrtoir.a v;Lt’i Fchocs 

7) Lov/-Order Corpenr*at3r Dcfjiqa 

3) Descent- P'aase I’ontrolloi" for 4 -d jation 

0) Irt init*; D’ aensicnal ?' ithe’^at ic.’.l Prcqrurr.’.iivj 
Fioblcns and ‘ inal Ci; atiol witli 

CciistrainLs 

10) Kobust '^on'^onsatot Oosien 

I 1) I.’ur.erical Mathcls for t):e I.y.ij ano/ iXpiat jons 

12) Poi tarlKition Methods in Linear Filt(.rinq and Control. 


/ 


1 
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r vvicr, y.FL 


IT’urin'} this tiro porioJ. the follo’/iufj reecivc-i partial salary support 
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Prof. N.R. ‘^andoll 
Dr. S.B. Cersi’.win 
Dr. K.-P. Dunn 


Students 
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F. I. IX 
S. Vareus 
R . Pueno 


P. lee K. Tonq 
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c . ir Tt". Arri ^^MS 

1*10 n.ijor intoracMon darin'-j fhis tino i;oriod bot’.scen M.I.T. {orsonaol 
and riASA/ATies ronoarch staff occarrcd dnriraj tho MA^A Ar.os Workshop 

on Syst*^n Fcliability Issues tor I’’.»tairo Aircraft, Tnis ri.'otinrj, ht.'ld at 
• f j-Qjr| 7\uqust 17—20, brou^lit totjothor a nuirlior of preninont rosoatv-hors, 

includin'! many from 'lA.SA/Amos and M.I.T. , to riiscuss a number of asjiocts of 
thoii' '.’ork that ir^pact on the question cf reliable aircraft; des:qn. Some of 
the v;oik outlined in Section D — specifically that related to reliability 
and failure detection — was reported at t'aat meeting. 
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D. '”?chni>~H, risc’ iifi ion 
1. P.~‘l ia>; i. M ty Opt i ni/, it. ion 

To dot'’, resoarch by !'r. Poujlas Hirdv;r»ll opd rioro.'.ior Miko Athans 
has bo-'n coticorned with th'’ ’.('proaontation of systo.-ns for us*; in I ’no studv' 
otT systons roliability and optimization. A tncor'itical rrar.owork is nor.fjht 
vihich is gonoral croufjh to ronsiJor mnv' Ji fforcnt tVT*^-5 of proLloins, V'-t 
is spocific cr.f’uqh to allow optimal d'’siqn of reli.iblo <>ystans. 'Hio 
frarrovork currently boin j ir.vo^-.ticjatc J accovnts t'«3r a] 1 pomihlo discroto 
failuros of a systom and allows f<^r optirir.at’ on ovar a s«:.t of system 
con*‘'iqurati'‘>ns v;hich ir.«?o,t a sr.'oci’^iod porfonnancn critori'-n. Given a 
sp^ecific problem, this .pproaah sooms to work well and has jivon rjood ri’sults 
for tho probLen of sensor allocation for a I'noar : ystom, wh'jro failure 
constitutes the lo.ss of '-ome set of sensors. The anproa'sV. also seems to bo 
api'l leablo to other areas, s.uch as the ch.cicc of linearizations for use 
in ♦’h.'a alaorii-i'ns for 'rvation ami control of siqnif icantly nonlinear 

systems. Some problems with the ap-procach include the thjo’ocical and 
cemi utational complexity encount'ore l in tl e application of the method. The 
theoretical complexity ray bo reduced ia the future with the addition of 
qeneral results applicable to larcje sets of problems. V.’a'/s may he foumi to 
redu<' > the computational complexity throuqh the careful use of approximations 
invo’vinq the set contai.ainq all possible modes of failure. 

A b»"apch of research which has b<’cn itiduc’d by the above work involves 
.1 study of multii'lo model adaptive cor. -rol of a nonlinear .sy.stem. A new 
ap! roach has been outlined which considers the set of linearizations of 
iht? system functic.a to be a finitc-s’-.ate space. Tlie problem beinq 
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i n V(!S f" ■ !.it O'i i '5 t.r; fi.)] lows: Tt 'i c*'r.L:.‘'>l ny'^t.i;ir> u_-ir'{ sfyir.j ■ 'in .u‘V»C 

i . a[ { roximt inq flu* sysf* n \;iLh .oro . ubstt of available li tioari-ial ioas 

ii at sonii timo t, liow sb.oulJ ‘"ho control lav; choo ’.o v/hat of 

1 i no irizatiojis S' Lo use at nooo future tine t'? One approach to tlui 
problen is to estimate U o "distance" betw.‘e:\ a syst.^n ami a qiven 
linoari.".aticr., aiul clvxir.a S' usinq <\ variation of the nearout-neiqhbor 
alfjorithn. W'.eu distance is dof'.neJ as the norm of the vector difference 
biitwoen the system state and a point, of 1 i near isatb.n/ t,Ti:j approach m^iy 
yiolfl a locally optimal method for ch.oosiiK; S'; however, th**re may le a 
I incar ittition which is a larqe distance .*^ron tlve state estimate, but which 
is quit.? close to the «har.i.-t >ri sties of the actual system. It would bt» 
nice i tboro existed a laasure of distance v.hich yielded a globally 
optimal choice S', t.'herlu'r or not this r.ea.j.ue t'xists i-. a question 
srill to be ans.-.ered. 


'I’l’is rese.irch 1; a specific c.ise of a more qcsieral area of prcbler.s. 
Stq'pose the researcher has a sv.sten v/ith a co;iti.v.:ous state space represon- 
tacion and .is interested in some coi.pact suijsot of the state space. If he 
considers a covoritiq for this subset as a finite-s^atc space, with links 
between stater: reprosontin'j the existence of a path from one elenirnt in the 
coverinq to a i-oint within amathor neiq’hborhcxji, in v.hat way docs optimiza- 
tion over this finite-state space corre.;pond to optimisation over the 
continuous space? Some connections with previous work appear in the 
discretization of a state space for diqital iurixjses. Future ap[)lic \tions 
ray inclr.de problems such a.s u.slnq this vi'-wpo.nt to qonerate ;ood, fast 
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11 ; r oximat it.ns to tho noliit ion oL o|>t mal control (-rol leir.; m l rcl iabi I Lty 
oi t i ni :: it Lon probl. r.L; conco’nocl with continMoun failuro rr loj in systems. 

2. Tiilnr Po tpi:*~ i on Algorithms 

Prof. A.S. willsky, Hr. S.n. Cnrshwin, Dr. K.-P. Diinn, 'Ar . K. Chow, 
an! Mr. R. Dnono have centinuoi tho study of fail iro dotocticn sytcin dcsi'jn 
methods, as outlined in the procedinq proqrc ss n i virt. Tn thin t ino i)ori.od. 
Pro-'. V.'ilLsky completed a survey^ of a numh>er of failu e dotoetii'n rothc Is, 
outllninq thfiLr 'trengths ar; 1 v.c iknos!,-’.; and the variens k*iy tr<i ’eoffs tn 
the designs. ’'n cdditloe, tiie 'II.P m.etho.i (soo the sur/oy pai er'' or the 
I rcvlous progress coi>ort) has 1 ^en studied in cro it ’etail. ^'oeei t icul 1 y , 
n extensiv; simul.ition st’idv» ’ ■ ' pg the siir.ili'iod tve-u imcnsicnal lonqi- 
tu.iinal d/n nnics of the r -dC airc* ift, hr; I'-.en ccngleted, and a r.*port 
devcril-ing tl'.o gual i' ..tivc he:’ avior of f.'.e d:teehor, filse al v: ’-coi rect 
detection tradeoff-., etc., is torthconinn. in *. his report wo al o •> icrihe 
a detailed anal’.tioal package that lias hcc’i anl is contitv.ing to he dev- loped. 
This package allows one to calculate probabilities of cor-ect detection, 
fiilse alarm, mis-classif Lcation of failures, and the exp-octed d.elay nine in 
dot cticn. 17e have afplied those techniques to the F-oCaircra’'t probl- .i, and 
have dote mined sc'/eral of the relevant tradeoffs between theso »erformance 
[larar.oters . Purt’nor development of these mcthoils will be don*; in th-i coring 
months. Specifically, the -llov/ing issues are to be considered in the 
irr.odiato future: 


1- .\.G. 'JLllskv',' "A SuTvey of Design M.;thcds for the Dot-'ction of Failures 
in Dyn.amic Systems," presented at the M.I.T.- :.7FA Ancs Workshop on System 
r.L'lial ility Issues for Future Aircraft," August 1^75; to be sol.nittod for 
publication. 





(1) The uot-. rnin.it ion of the n«'r.jitivi ty of the dntcJtcii vii ::’nulation'» 
of tlie F-Hc aircraft vith orrors in the rel> v:»rit aoro coeft’ if-ienti*. 

(2) Development of <i ain.plifi’ d detection syaten and .1 Mtudy of its 
per ♦’'ormiiiCe via analyijis and sinulations. 

(3) Cl- velopr.ont of detection retho<lu for other types of * ai lure modes 
th m those presently beinq explored. 


3 . An alysis of Koni i near Filt-^r n 

The work of Prof. «\.3. Will'ky arul r!r. r.T. rtarois reported in the 
pieviou.s research rei'ort has been wixtten up in nover.al p-a{>or-> (sec the 
references at the end of this rej.art) . V.’or.< in the .trea of (.online. ir 

stochastic fil rs has cc.ntinuec alonq the follt vinq linos: 

1. V’o are {ye-.s ntly ir.v-?stiqit inq the ex ton. Lon of the cptiruil nonli.iear 
e.sti-.uiti-f. results o.' ft’.rcus and '..'Ll!sk/ {.;er> tl;e 1 1 : ference s) to 
obta*..i .accurate s’ ’ optir-*! estirr »tijrs ''or other classes of nonlinear 


sys terms . 


•y 

2 . In his thesis', riarcus dev.;l<-ped .1 r.ubopl inal cr.tiration technique b.ased 
on harmonic .lualys’s. This techni p:e is b 'inq applied to desiqn 


estii ators for attitude e-’tim.ition system;. 

3. Prof. Willsky .and Vr. V. Klebanoff are studying tl.e "r.armonic analysis 
filt. rs” d .veloped by VJillsky an.l J’jircus in order to d itermino 
analytical n.-asures of their perforr.»nces. 

4. Prof, h'illsky and Xr. J. Ktorno are considerinq the problem of the 
optimal demcdulition of a .si .'.usoidal siqr.al. The re the 1 here is 
quite different from tl.e Fourier series methods of .’’ircus and h’illsky. 



Thi:; ru- ; »i ji-Oiioh involvor. a : i ■ ♦ co:. LloiaMon of t!,»‘ so-r.ill<l 

” r*-M'o.'i«‘ntat ion theor**'T” i-xi lori f>'i tiiu cor.;I it. tonal prr>l*abi 1 ity 
.lonsity. We an* proRoiitly i'e»-ivin j a; proxiinat«’ ex» i for the 
exj T>ctat i ons to evtiluate t he cttn-iity. Thin at pfo«iCh r? -.y prove 

to Le .'x^.rrrely I’sefnl in a wide variety of nonlinear I’ntimat ion 
ai'plicjitions. 

5 . t.’ninq tha sare Volteria sei-i-.j n.»thcd t an Mioie ev.ploril by Marcus^ 
Prof. 'I’illsky in sti*'lyin<i tho rtaljil’ty of lit. ear t>- -Ihack !yit ns 
cont’.iiiiiej randomly varyln*} jains. 

A, o f r • T i r 'I Ti' e f'etays 

Prof. Uillnky and *Mr. K. ;'’..o;;7 h tvo ar.dirtakoji a study of the 
of f-'otibaek oor.. >cn.sators th.it i.ict' 70r.it? tire d lays. <.'rinsidor the 
linear systen 

x(t) ^ «'j<(t) I- r'i(t) (1) 

y(t) = Cx(t) (2) 

u1;ore ;<(0) i; arl itraiy. Tf w.? t. strict ourselv-s to linear facdb.ick laws 

u(t) “ (3) 

the Lost wo can liopo for is an exp jnential decay of the output. However, 
if we allow time d'l.ays in the control 

u(t) ^ K^:<(t) + K^xCt-h) (1) 

one oan iften dri .e the output to 7 ero ( .nd k‘»ep it there) in finite tine. 

In his thesis^, ."r. Kwer. j has consi ’ .'f. i a develoiner.t of such a desi gn in 

2 . . . 1 . ;!aicu. , ''list iruat ion and An.ilysi.s of t'onlir.far rtocii.astic F'/.-jtem.s , 
rh.D. Thesis, Oopt. of Electrical F.noinoering md Computer Science, 

M.iy 107 S. 
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l-.,i sin.jl.-outrut cai:o c.uvl has (h.viseh i iLOceh.ti for •U-rr.i n iivj Y.^ ml 

w' ' n (A,^) IS a controllable pair of dircnsion at least 3. Or - of the Laaic 

4 

r '-ult.'i of this -Jttidy is ref;orti.d in a brief :;umrary laper. 

The next 5 t«’p;» in this -.tuily are: 

(1) The considcrwMcn of the dual problem of •'leailbcat observers" — v.e. 
asr ur.in.j (A.C) is obac.vahle, desi.fn a syr.t n of the form 

^(t) = /v:^.t) t nu(t) 4 K^lv(t)-cx(t)) ^ K 2 (y(t-h) - c:(t-h)) (3) 

j.i that t he e:: -irition or-or 

u(t) - x(»-) - x(t) 
i i .iriv'en to ia finit^^ 

(.1) I xt. r.a Kwon j’.'. r •- ults to the nultiple cutp'it t:‘* 

(}) T.- 'lclei .h.l.iys i.'. f-.o sv'-t . -Ivi -^ics a:vJ cuH-iit. 

( \) Con i*h r U.e d..si ;a ‘ t co . p. nsato. J of the forn (5) pl-ier. 

ult) - .N^x(t) + X^x(t-h) 

JO that systcra output is null<-5 ia fl.'.ite tirvs. 

^;inca r..ny jhyaic..l ayat. :aa hav inhercat tiro delays, the above study 
jaouU aid ia our oveiall undoratandii. . cf these sy. tcns ai.J in our ability 
to vleviso effective control cysters for them. 


YTITl! . - ' ^^(eeraic btT icture an i F.t;rotioa of lire Delay L/stems," 
rh.r. Thusi-, Dept, of rloctrical Kn^ii -eriny ani Coirputer ncionco, 
C.ai'hridae, I^tiss., r«;:p't. 1375. 

4. E.-.-S. Kvom ,.r.4 A.s. •■Flnit»-vir.c :<=ro Error Co.-rol Vi4 Dolay 

Focll-acl:," U75 :ni;E Cont. on D-.-clsion ,vid Control, Ecistor., Doc. V>^S. 
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^ . r i T t t.n? < "<;»-! < (.f . r rv»«:* in 

Prof. Willsky .^n*l Mr. M.K. P\>b*.rts ht>vo con? iautHi offcrto in 

rtLton’|’tin<’j t»> imiler'? nnd th»i fun<:an.'iU.»il 1 •.;ritaficn i of inploru'n- 

tat ion of ironLrol uy-rtems ( '.ce the pr^.<'<*din»i f roqrcsr rej'ort) . Tlio Li^’ic 
r^ntivatiosi for this study is to take a step in d*n'-lrpin»j a control 

dc .itjn m*3tho-lolc«jy that directly incorporates cor.; utt-r t r dvlocffs. fticcifi- 
cally, \.’o h.^vi» b<>cn con-: idc ring the design of digital controllers for linear 
svst'.TS based on quadratic p- ’-fonnanco no.i:;urf'S , C<ir study has incbhlc.d 
OX! licitly f’l** constraint of a finite ti.no to p'}r.''/tn a ’•c'l alf.ipl lea? Icn, 
and •.’o have concentrated rn finding ’cw-orlor exarples for vliieh ainpler, 
fa.?ti r, "sul opt inui I " control lav.*s ou^per the slower " [ ti^il" coitrol, 
whi'h roguiros a nuch slower sani>ling rat*. Thr* i v-nt.^il aisi of tain stuly 
i ' to <lev..lop conco['tJ ^'or th * centr of large--. cile s>- ti-r':, nd onr 


initial ideas v/i’l be described in detail in Xr. 


' ikorto* forth* oriing th*.-iis. 


b. r >tir».>t ion for '-v- *-er.- Kchoos 

A problen of gie.nt practical ir.p-or ta.ice i'* the esti :'a-.lon of signals 
which tiriy he dolny-.J i.i time or which r*ay 1 •» received to-;-^ther \;Lth delayed 
v-^rsions. T'ao siirplest such r«i<lol is 

y(t) « h^(x(t)) r h^ixlt-T)) + N(t) (1) 

vlicre X is t!ci signal to be recovered, M is receiver nois- , and T is a ti^'e 
delay — cithe»- known or unkno--in. Prof. V.’illsky and Mr. P. Kwong have studied 
tiio ostinatirn of x given such a received signal (actuilly a n»ore g*'neral one 
involving the possibility of distribetc-l delays) when the delay is known, i:»d 
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r,o'-‘-o !»'•• ' :i >t.or D»;n i ()n 

Tho offorts oC Trof. T.T.. Johr.fjon , Mr. K. Tong, Mr. li. ’'..lorenrj, .\n\l 
M.r. li cont.inuG to locus on the dusiqn oL' lov;-ordor dynmic coirpen- 

s.itors for r.ult ivar iablo ILnoar systems. The interplay between forruilat ion 
of the o; tii.uil f'ontrol prohlctn, existence of fiolutlons, and structure of 
t*'.e optinaL cont»*ol law is tho focus of thi.s research. 

In a fort'.-.coninq lldt Technical Meir.orandun v/e ’..ill pre.sert tho solution 
i f the roduced-o-'.ler r,ljt,i>rver-hTicd optimal conpensator design problem for 
I’.e st.u'.danl rua.lratic perforpance index. Reduced-order observer-based 
C' g ens >tor.s (’.o., tho.io based on a Luenber«j2r observer './tilcin estimate.s 


vh'^'.e ] i ;.i-ar corbin itions of states roqv.lrcd for the optimal input signals) 

■f. tve Lein stuli I by hrasch and Pearson, rortniann and Lilliimson, F.oi-hsch i]d, 
•1 .‘Lira -hik, ocl’.’'‘rs. The mi-iunal dim.ensicn of sudi compens.ators Is 

' r. ’.n to I: t’;.e : u<im.u:i of ( i) th .2 cl ■. v.djilit/ Li.dex of (A,C) minus one, 

• I (L) (n-.n) r/m ••he re ’ nurb ’r of input.i a.nd n== order of pl.mt, and m= 

r; be»- of o’.p j ’.M s . to-. -r order linear cor.penritors '.vhich stabi lise the olaiit 

” iv .tin exi t, Li’t tho’/ caniiot be b.\s«.d on oL ior’v’cr theor'/. The above 

• o’ or.; pr.i- .il!’/ use \ pol-i-pl .iceinent approach as erpoased to ninim.ir.inq 

t!..* ;.ta ’-.rd ::.iL>-atic i ■-*r ‘’orr. ncc. inlex. We shew that tho oi '-imal structure 

I . ’ r.d • I ;ivc n ly .\n t '^r’/-‘r ’.Lich estimates (Fx) vd’.ere F are the Itilman 

-ii.ns or th'. 'ii-.ear i .ulator I'r.osler', and qi'/e an explxcit prcci^dure fer 

p n ‘ru ’flo.i of * he rp* .'.r* ; ;ns ito»’ . i^nalcso-ij results were proved b’/ 

6 7 

” 1’- r i • d) , i i!l lio in (1073) f.^r tl j (n--m)th orde-r observer. 

W r v t .* 1 . tc ’'.'I* 1 • 1 • .o’/id'! d«;.;i ;n nt^theds for any order d.j'*n to 

• a.-- r : — :m c * ’ *- 1 aho’ e . 



V’o havo uaod init.iai-.'jtat'.o avora'/in<j over .i rair’-on initial plant 
t:o detamLno tJ'.o opMrUii initial statu of .such coi puusator.s . This muthcd is 
not entiri!ly sati sfactoty , and the oiitinal coinpensator clynanics in fact 
depend on the initial r:tate covariance in tlvis case, A lar'/o covari '.nee 
iir, plies faster optir.al d.yn.anic.s. Tiicrc; appears to h-o an alternative, an«1 
scraev.’hat more appoalii\c; irtorprotation of our necessary conlitions; an 
oq>iivalunce relation hotvoon the randcin init.ial state and th.o stochastic 
steady .st ite p<roblcrr.s can be e.stahlished for t hiU (n-m) order corp'-nsator. 

In the stochastic steady st.ite problei.i v.e a.iSUi.'.G t’.’.e same plant and coirpen- 
sator structure, except that the filaut is driven by a randon (white noi.se) 
input and i.T in stcch >,stic sLua ly .state: \.e seek to minimize the expectation 

of the inteqr.u'd of the usual <iUad.ratic j-er fr/S.v.aiiCe index, which i.s tirr.o- 
ir.v.iriant. Perfect observ-itioirs are -i.i'c.iiie .1 . l-^sor. rally, the rol r of the 
Initi.rl state covari .n.co is ; bayed 1;/ the injcit noise intcn.sity of t'.'.is now 
problem. 

The pro.bb.n of sens' ti’/ii.y analysis and inpjlerr.v.itation of such djsiqns 
are still under invest' i rjatio^' . V/e inteml to develop .a rod.uccu-ordur desiqn 
for a realistic aircraft r'.cd.-'l (see L :low) . 

Some authors have clained th.at ob.sorver-basoc! dcsiqns a»e s nsitive to 
r.rodcllinq error and al.;o Ja poorly in req'tlatinq stet.d.y- state errors in the 
face ot disturba.sces. They proi ascd ir.cl'idinq fon/ard-loop inteqrat ion in a 
rather ad-hoc manner which mimics convcntioucal irequcncy-domain methods. We 
qrant that model sensitivity may be a pr.,blom, but expect that (in the absence 

6. Miller, tb.i., ''S[e<;i'ic Ci tinai Control of the L near Pequl.ator u.iinq a 
r'inlm.tl-Or Icr Observer," Int. J. Control, 1073, Vol . 13, No. I, pp. 130-150, 

7. Pi invillain, "Optimal Compensator Structure .'or Lincmr Time-Invariant 

Plant '..'ith InacciV.sible State.''", bSb Report , May 107*1. 
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of r.io<2cl niirrna-.oh) tho structure of the 1 Lnitol-cUmcnsicci conpensetor shoul,<i 
l)L- .sufficiently fjonoril to allow for fc^^Viird loop iutoijrution to inchiovo 
zero steady state error, when recpiLred by tne quadritit; perfornidiicc index. 

This deironstrat ion retpii res a sliqht goneralizatron of the problem Lomulation 
above to allow exogmious inputs to the compons.itcr . V?e are considering such 
problems . 

Kc also continue to investigate the (luestion of uniqueness, or the free 
d.osisjn pararotcrs of the cptin.»l compens>\tor . 


n. p.^sc' nt-'''has • ; Control l er ' 4-0 M.r/i-gation 

In a recent thesis bv P. Lax wo have investigated the feasibility of 
d'':', igning a t iT.e--v.iry iisj control law to aahi*.*ve simul taaoous Irajectoi/ follov/ing 
and .iccurate control of r iuway -i.'.ark.er arrival tire for the descent ph.^so of 
fli'^jht. The procad.jre I'or 1 ir.ec riti a j t'.aa r.o;.l.inear alroratt cquat lot\s of 
motion is v;e.il -known, but is seb'.oin carried out in a way con.; i stoi.t with the 
"L9C; .losi-n philosophy". In addition, the descent phase cf fU<iht is usually 
liundlod by gain sc';eduling models, ratli*-*r than '..esipning true ti.r.e— .’ar/ing 
optimal gains as csi'>«au.;ed by "‘"he" theory. V.'e h.ave carried cut all these 
cliiugs in a consistent manner for the lcngitudi;'.al d.ynanics of the Boeing HJ 
on an a{ pro.ich to Ixigan Aimort. Tl.e time-varying liucarLzcd controller was 
furthornore implenirn! ^'d on the r.onlii’.ear digital sinul.itii>n ana pilot-tested 
under various mean-wrad conditions cn tne ESL cockpit simulator. ihe result.! 
were quite encci raging a.nd will Le reported in a forthccning I.£L Report. 


0. Infinite Pit - nr. ion/ 1 Mather P. ical r >~r a ram’- i r a rnblems atid Optimal 
Control ?>-rd.lo:n with Con.itr.', ints 

There recently has been considerable progress in non-linear progr.i ruling, 
both theoretical and conputatien.il. It is recognized that duality plays a 
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C"!ntr il rolo .u;il du ll TnGtlio<lf3 t‘ » y turn out I b»; tlio r.ont ip.|'-ort.\nt comi-uta- 
ticnally. The uuo of nuiinontud Lrpanqians (1) 13 cxi <?ctc;d to play a kuy 

rolo. 

Opon loo{) optur.al control probloirs are Lnf inito-diri*insicnal mathonatical 
proqramning problcns. Howovor, they have conuidcrablo structure (niunoly, the 
time structure) v;aich i.’ust bo exploited. V'ie'.inq then as infinite dimensional 
mathcnatical proqramjninq problcns has the advantage one can focus on issuofj 
which n:qht nou be evid*Mit otb^rwise. It has lonq lieen felt t’aat a synthesis 
of mstln.m.atical prograi.r'.inq and optimal control would bo aa inportant devol- 
oj'nent. Kurthermoro , it has also been felt tliat the so-called dir«>ct notheds 
(Rayloigh-Rif;) , whioh have been largely nc‘ilcct«' ! in conputatior.al wor?< for 
cptimal >'o:it‘ol problems,, should allov; us to use nather.atical [rcqramr.ing 
in an efficient way (t'lat is, exploiting the time structure). 

In recent i:or). {[2]- [5]), '.,e fed v/e have taken important stops in 
tliis directicn. .‘Vjro { reciselv, we have given a rc asonably comp lot,* I’uality 
theory far convex c[itiir.:.l control problems and also shown h.o*./ the flual problem 
can bt; efficiently solved using the Rit:; method. 

u’ork is II' iw progressing in several directions: 

a) L's I’.'.g similar i leas we arc investigating the following linear problem: 

M’ nini 2 :e 

1 

c(x,u) = J (a'x(t) -t 3'u(t))dt 

o 

rubject to 

X - Ait)x(t) + D(t)u(t) 

K^u(t) < 0 f K^x(t) ^ 0. 
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V/.J cotijc^-tufi' that, wltv.out. a '.li.xtQr t./i.e condi.t.ion ;mcl wlt-.li fc isibi Lity 
or: rhe pri .lal and dual problnmj that-, trhoro is no dualiry qap. Ihis i ivnilt 
will ('onoLMllno an oarlior result of Levinson [7] and ';ill have qrcat conpu- 
taticnal yiqnif IcancG. rrc.-blemn of this typo arisf', for enamplo, in Dyitamc 
TA.‘cntief fiod.cls and in certain prohlt^ns of cor.nunication networks, 
b) After d i scrot i:iati on of the dual proLlt-'M uuinq finite ei'_ .lont motho<l:> 
we >re left with a T.athcu.watical prograrininq T'toblcn. In fie quadratic 
c..:;o we are left v/ith a quadratic prof/rai ir.inq prohden. Tt ai i-o irs however 
that those (niadratic proqraraminq i roblcns chc.*ld not be solved using 
coi, '■‘lerv'ntarity faeoty. lie arc investigating the .structure of the m.iMiona- 
tical I rcgrarraii g probl.- with a vie'./ to obtaining efficient tec’.miquen 
((•orh.ipr. thicorpon iticn r’.i.'. thods) for th.t;ir solution. 
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in Conve:: D^'o-jr i. .'rine, " .ST7. ! J. <in Concr .l, pp.2Cd-2flj. 
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10. C'- H cnGntor 


Mr. P.K. '..’01U7, Mr. M. ?!afanr«', and I-ro' . M. Atl-ans hava boon i nvost Iqa cincj 
s-'v^ral anf'oct:; for dool'in of -Ounlrablo Conbiol .a dositja and 

Kalr'.an . ilPcr <i*isiqns that ara t undor large { ara:ioter variations. 

In the dor.inn of comf.onsat(,i s for the control of dynanical system 55 , 
it is ccrrjnon cnjineoring practice to introduce the asnunption that the 
dynamical syston to be controlled can be adequately n.cdoled by a finite set 
of ordinary linear differential eiuations. t'nder itiia aesunption, t.he 
engineer can, with r^^lative ease, synther.ir.o a linear state*- CoedLack con- 
ponr.ator for the model which will r'inirri^.o a qtiaclratic perfermance ii.dox. 

The class of cempennatorn pcc luoed by this design proco inre hav^i cone to 
l;o g'inerally I'-novn as linear-qua'lratic- (bO) controllers. 'Tne research 
carri-d out is aimed at d-t jL'-.ining criteria ly wliich it will be i osr ible 
to judge a orio» i whether this design prodecure can be ertpec'. ;d t'"' yielfi 
g'^od ro:;ult3 for any i.ittica’ar d/numcal system. 

F'evilts reported in the lit'‘raturo to c'ate on ta.e validity of this 
design procedure have been largely limited to a p»osteriori tests of the 
soa.sitiv'ity of the cf^ntrolloil system to iufinitosimeal per turbat -ons in the 
’•aricus parameters lioi-tribing the systsm. However, very little is ..no\/n 
rp’.antitatively about the validity of the use of LQ controllers. One notable 
exception is that ^ quant Ltativo hour. !s hnve been found for the 
ninirum quin and. [dnase margin of single- input/singlc-output lirioar systems 
with I.Q controllers. Recently, the gain margin (but not the phase margir.) 
bound for L<? controllers has hem generalized by P.K. Wong to include 
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rrulLi-loof. linear feedback systems. 

We are presently cxair.inincj analytical techniques which are expected to 
make possible more widely applicable quantitate statements about the validity 
of the use of LQ controllers. Currently, efforts are beinq directed to the 
generalization of the aforementione notations of a priori qain and phase margin 
to LQ controllers applied to nonlinear systems with multiple feedback loops. 

The immediate objectives are to 

(1) provide a theoretical foundation for the empirical fact 

that LQ controllers tend to be insensitive to modeling error, 

(2) provide engineers with moans for determining quantitatively 
the magnitude of ino<leling error which can bo tolerated 
before it is necessary to question the reliability of a 
system with an LQ controller, and 

(3) permit better understanding of how such effects as phase 
shift, time-delay, and saturation can bo expected to 
affect tho performance of multi-loop LQ control systems. 

ri*el iminar>/ results indicate that a quantitate bound on tho amount of 
nonlinearity and phase shift which can bo tolerated in the feedback loops 
of systems with LQ controllers can bo determined a priori for entire class 
of systems with LQ controllers. 

Tlie results for robust LQ controllers can also provide valuable insight, 
namely the design of Kalman filters for which the covariance of the fake 
white noise can play tho role of a design variable. In the long run wo 
ho{’o to obtain a unified quantitative design procedure for the overall 
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LQG robust dosiqn concept. 

Faculty M. Ath.ins 

Research Assistants 

P.K. Wonq 
M.G. Safonov 

Reference ; P.K. Wong "On the Interaction Structure of Feedback Control 
Systens”, S.M. Thesis, K.I.T., August 1975; also reput ESL-R-625, 

September 1975. 

11 . ?iumerical f*ethods for the I.y»T;"mov Equation 

The importance of the algebraic Lyapunov equation in problems of 
steady state estimation and control of time- invariant linear systems is 
well-known. Mr. T. Athay and Prof. N.R. Sandell Jr. hav’o initiated a 
research effort aimed at development and careful numerical analysis of 
algorithms for Lyapunov equation solution. The research will be docum.onted 
in the forthcoming Engineer’s Thesis of T. Athay (11. The following 
preliminary results have been obtained. 

The notion of conditioning is basic to numerical analysis. V7e have 
obtained a precise notion of conditioning for the Lyapunov equation. We 
have been able to dem.onstratc that the often heard statement that systems 
with widely spread eigenvalues have ill-conditioned Lyap-'unov equations is 
correct in a precise sense, and moreover is a property of the equation 
and not the solution technique. 

Although different algorithms for Lyaj unov equation solution cannot 
affect the problem's inherent conditioning, they certainly will have different 
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error 5 )ro! oi ties nn>l operation counts. V.’e have analy?.ed a nurber of the 
more popular alqorithms with respect to these properties. 

Finally, motivated by the bloch coupled structures conmonly found in 
larqe scale enqineerinq systems, we have developed certain alqorithms that 
explicitly exj loit these structure. ,\nalysis and numerical test of these 
alqorithn is presently underway. 

Peferences ; 

1. T. Athay , Numerical Analysis of the Lyapunov E^mation with Application 
to Interconnected Power Systems. 

12 . Perturbation Methods in Linear Filtering and Control 

Large scale aerospace and other engineering systems are character ized 
by interactions of smaller subsystems [ll . These interactions may be 
weak so that an approximate analysis of the decoupled subsystems may bo 
employed, or they may be so strong that certain subsystems can be aggregated 
into an otjuivalent subsystem with fewer state variables. The mathematical 
tools of perturbation theory are an appropriate means of analysis for such 
situations. 

Mr. D. Teneketzis and Prof. N.R. Sandell Jr. have initiated research 
aimed at extending the existing tl'.eory of perturbation methods for deter- 
ministic optimal control problems to stochastic filtering and control 
problems. The research will be documented in the fortlicoming Master's 
Thesis of D. Teneketzis [2]. The following results have been obtained: 

For the weak coupling problem, the analysis proceeds via regular (non- 
singular) perturbations. Conditions have been obtained for the stability 
of filters designed for the perturbed (decoupled) subsystems but applied to 
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the actual, couplod systom and conputable bounds have been obtained for the 
resultinq deqradation in performance. Analoqous results have been obtained 
for the stochastic control problem. 

For the stronq coupling problem, the analysis proceeds via sinqular 
I>erturbation theory. An interesting and nonintuitive result characterizing 
the behavior of the optimal filter as the coupling becomes infinitely strong 
has been obtained, and the analysis has been extended to systems with three 
levels of coupling. 

[1] N.R. Fandell Jr., "Information Flow in Decentralized Systems," 
Conference on Directions in Docentral izrd Control, Manv-Pcrson 
Optimization and large-scale Systems , Sept. 1975, Cambridge, Mass. 

(2] D. Teneketzis, Pertiirbation Methods in Stochastic Control , M.S. Thesis, 
forthcoming. 
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E. F I\’A:.'CIAL status 

As of December 1, 1975, .inproximately 301 of the qrant funds ($70,000) 
wore exi anded corresponding to 25% of the qrant period. 

No consideration for cost sharing is included in the above figures. 
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